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ABSTRACT 

(N ; 

^ Context. Precise S abundances are important in the study of the early chemical evolution of the Galaxy. In particular the site of the formation 
remains uncertain because, at low metallicity, the trend of this ff-element versus [Fe/H] remains unclear. Moreover, although sulfur is not 
bound significantly in dust grains in the ISM, it seems to behave differently in DLAs and old metal-poor stars. 

Aims. We attempt a precise measurement of the S abundance in a sample of extremely metal-poor stars observed with the ESO VLT equipped 
with UVES, taking into account NLTE and 3D effects. 
Methods. The NLTE profiles of the lines of multiplet 1 of S I were computed with a version of the program MULTI, including opacity sources 
from ATLAS9 and based on a new model atom for S. These profiles were fitted to the observed spectra. 

Results. We find that sulfur in EMP stars behaves like the other a-elements, with [S/Fe] remaining approximately constant below [Fe/H]=-3. 
. . \ However, [S/Mg] seems to decrease slightly with increasing [Mg/H]. The overall abundance patterns of O, Na, Mg, Al, S, and K are most 
^ . closely matched by the SN model yields by Heger & Woosley. The [S/Zn] ratio in EMP stars is solar, as also found in DLAs. We derive 
' an upper limit to the sulfur abundance [S/Fe] < +0.5 for the ultra metal-poor star CS 22949-037. This, along with a previously reported 
rS measurement of zinc, argues against the conjecture that the light-element abundance pattern of this star (and by analogy, the hyper iron-poor 
. stars HE 0107-5240 and HE 1327-2326) would be due to dust depletion. 
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1. Introduction 

The determination of the abundance of sulfur in the atmosphere 
of old Galactic stars is of interest for two primary reasons: 
-First, since unlike most other heavy elements, sulfur is gener- 
ally considered not to be significantly bound in interstellar dust 
grains, it is a good indicator of the chemical composition of the 



* Based on observations obtained with the ESO Very Large 
Telescope at Paranal Observatory, Chile (Large Programme "First 
Stars", ID 165,N-0276, PL: R. Cayrel. 



ISM and thus allows a direct comparison between the abun- 
dances in the DLAs (damped Lya clouds) and early Galactic 
matter 

-Secondly, the formation of sulfur remains controversial. Sulfur 
is generally considered as an a-capture element such as mag- 
nesium and calcium, formed preferentially in massive type II 
supernovae, including those of the first generations of stars. In 
analogy with other a-elements, [S/Fe] is expected to be con- 
stant and positive at low metallicity (see e.g. Frangois, 119871 
I1988I I. However, the trend of [S/Fe] versus [Fe/H] in the early 
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Galaxy is still debated: Israelian & Rebolo (1200 11 1 and Takada- 
Hidai et al. (I2002I I found [S/Fe] to rise with decreasing [Fe/H], 
while Ryde & Lambert (|2004] | and Nissen et al. (2004, 2007ll 
found [S/Fe] to remain flat. Finally, Caff^au et al. (2005, 2010a), 
combining their own data with a large sample from the litera- 
ture covering the multiplets 1, 3, 6, and 8 of S I(see Table [U, 
suggest a bimodal distribution of [S/Fe] for [Fe/H] < -1.0 to 
explain the high value of [S/Fe] found in several metal-poor 
stars. 

Table 1. Parameters of the sulfur lines 



with 5 pixels per spectral resolution element. The spectra were 
reduced using the UVES context within MIDAS (Ballester et 
al. l2000l l. The signal-to-noise ratio (S/N) per pixel in the region 
of the S I lines is typically ~ 180 for the giants and ~ 90 for the 
dwarfs (or 400 and 200 per resolution element, see Cayrel et 
al.. 12004b . 

Three Ba-poor giants (McWilliam et al. 1 19951 1 were also 
observed as part of the programme, under similar conditions as 
the other metal-poor stars (see also Andrievsky et al. 201 li to 
be submitted to A&A). The S abundance could be determined 
for one of them (CS 22949-048), which has been added to the 
sample. 



Mult. 


wavelength 
(nm) air 
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log Sf 


Xex 
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-0.51 
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869.4626 




'P3 


0.08 


7.870 
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675.6851 


4p-5d 


5P3 -^D? 


-1.76 


7.870 


8 


675.7007 




5P3 -5D° 


-0.90 


7.870 


8 


675.7171 






-0.31 


7.870 


IF 


1082.1176 




3P-'D 


-8.62 


0.000 



A rise in [S/Fe] with decreasing [Fe/H] is difficult to ex- 
plain. It has been suggested that SNe with high explosion ener- 
gies (hypernovae) could have made a substantial contribution to 
the nucleosynthesis of elements in the early Galaxy (Israelian 
& Rebolo 120011 Nakamura et al. I2001I I. or even that they 
may be responsible for a delayed deposition of the supernova- 
synthesized products into the interstellar medium (Ramaty et 
al. 12000b . 

On the other hand, Takeda et al. (12005b noted that, even tak- 
ing into account the non-LTE effects on the formation of the S 
lines, the resulting abundance of sulfur depended on the mul- 
tiplet used for this determination. There was, in particular, an 
apparent discrepancy between the abundance trends deduced 
from multiplets 6 and 1 . 

Taking advantage of a new model atom of sulfur (Korotin 
120081 12009b , based on new radiative photoionization rates, we 
here report S abundances for 33 very metal-poor stars, includ- 
ing 21 extremely metal-poor (EMP) stars with [Fe/H] < -2.9. 
Our study is based on a precise high-resolution NLTE analy- 
sis of the lines of the multiplet (hereafter Mult.) 1 of S I at AA 
921.29, 922.81, and 923.75 nm (see Table[B. 



2. Observations and reduction 

All stars were observed as part of the ESO large programme 
"First Stars", using the high-resolution spectrograph UVES 
(Dekker et al. l2000b at the ESO-VLT. The spectra were dis- 
cussed in detail by Cayrel et al. (12004b and Bonifacio et al. 
(12007112009b . The resolving power of the spectra is R ^ 43 000 



3. Choice of sulfur lines and determination of 
fundamental atmospheric parameters 

In cool stars, five multiplets of S lean be observed in the red and 
near-infrared spectral region (Table[T]): Mult. 1 at about 922 nm. 
Mult. 3 at 1045 nm. Mult. 6 at 869 nm, and Mult. 8 at 675 nm 
(Caffau et al. l2005ll20T0ab . and Mult. IF, the forbidden fine at 
1082.1 nm (Rvde. l2006l Caffau & Ludwig, |2007b. However, 
in EMP stars, the lines of multiplets 6 and 8 are too weak to 
permit a precise analysis (see Fig. [1]). The lines of Mult. 3 and 
Mult. IF, are non-detectable in stars with [Fe/H] < -2.9, and 
are in any case outside the range of our spectra. The S abun- 
dances reported in this paper are therefore derived from the 
lines of Mult. 1 at 921.286, 922.809, and 923.754 nm. 

The line at 922.809 nm is located in the wing of the Paschen 
( line, and the absorption from the hydrogen line has been taken 
into account in the computations (see Fig. |2]). This spectral 
region is also crowded with telluric absorption lines of H2O. 
Since no hot star was included among the programme stars, we 
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Fig.l. Spectra of BD-H17 3248 and HD2796, showing the 
fine at 869.46 nm in Mult. 6 of SI. We detect the line in 
BD+173248 ([Fe/H]=-2.1), but not in HD2796 ([Fe/H]=- 
2.5). A synthetic spectrum (full line) corresponding to the 
equivalent width of 0.26 pm found by Israelian & Rebolo for 
this line has been superimposed on our spectrum of HD 2796. 
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Table 2. Equivalent widths (pm) of the S I lines of multiplet 1 . 
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Fig. 2. Spectrum of the giant star CS 29491-053 ([Fe/H]=-3.0) 
in the region of the S Ilines of multiplet 1 (wavelength in nm). 
The synthetic spectrum (red line) has been computed for the 
adopted abundance (Table [3]), and for an abundance 2.5 times 
smaller (thin blue line). The S line at 922.8 nm lies in the wing 
of the hydrogen line P^; all other lines in this region are of 
telluric origin. In our analysis, all S lines blended with telluric 
lines have been rejected. 



did not try to remove these telluric lines, but carefully com- 
puted their position relative to the stellar lines using only S I 
lines free of contamination. At least one or two lines in each 
spectrum in general are free of blending. 

Table |2] provides equivalent widths of the S I lines, since 
they were published neither in Cayrel et al. (120041 ) nor 
Bonifacio et al. (2009). In this spectral region, the observed 
spectra are affected by fringing, which is not completely com- 
pensated when dividing by a flat field (Nissen et al. 120071 ). 
The continuum level is therefore often not well defined, and 
we estimate that the error in the equivalent width can reach 
0.4pm for giants and 0.5 pm for dwarfs. Since in all the cases, 
the sulfur lines in Table |2] are "weak", the uncertainties, as a 
first approximation, are proportional to the EWs of the lines: 
A{log ab) - A{log EW). However, these equivalent widths were 
not used in the computations since we directly compared syn- 
thetic and observed spectra; they are given here only for com- 
parison. 

In total, we could thus determine S abundances in 26 giants 
and 7 turnoff" stars from our original sample. 



Wavelength (nm) 


921.29 


922.81 


923.75 


;r=x(eV) 


6.525 


6.525 


6.525 


log gf 


0.42 


0.26 


0.04 




EW (pm) 


EW (pm) 


EW (pm) 


T\imofr stars 








BS 16968-061 


1.68 


- 


- 


BS 17570-063 


1.30 


- 


- 


CS 22965-054 


1.30 


- 


- 


CS 29499-060 


1.80 


- 


- 


CS 29506-007 


1.98 


- 


- 


CS 29506-090 


1.63 


- 


- 


CS 30301-024 


1.95 


- 


- 


Giants 








HD 2796 


4.80 


3.20 


2.76 


HD 122563 


- 


1.92 


1.29 


HD 186478 


- 


2.70 


1.68 


BD+17 3248 


- 


4.88 




BD-18 5550 


1.65 


- 


- 


BS 17569-049 


- 


1.20 


1.20 


CS 22169-035 


- 


1.00 


(0.76) t 


CS 22186-025 


(2.09) 


1.52 


0.79 


CS 22189-009 


0.72 


- 


0.32 


CS 22873-055 


2.37 


- 


- 


CS 22873-166 


2.00 


1.24 


1.08 


CS 22878-101 


- 


- 


0.72 


CS 22891-209 


1.28 


(1.00) 


- 


CS 22892-052 


1.18 


- 


- 


CS 22896-154 


1.98 




0.79 


CS 22897-008 


0.90 






CS 22948-066 




1.20 




CS 22953-003 


1.43 




0.74 


CS 22956-050 


1.08 






CS 22966-057 


2.49 


1.04 


1.15 


CS 29491-053 


1.53 


1.02 


0.88 


CS 29502-042 


0.83 






CS 29518-051 


2.30 


1.34 




CS 30325-094 


1.33 


0.98 




CS 31082-001 


2.19 


1.90 


0.97 


CS 22949-048 




0.55 


0.45 



t Equivalent widths in parenthesis are very uncertain 



4. NLTE calculations and 3D correction 

The atmospheric parameters of the stars (Teff, \ogg, [Fe/H]) 
were adopted from Cayrel et al. ( l2004l l. Hill et al. ( l2U02l l. 
and Bonifacio et al. ilQQli and for the Ba-poor star from 
Andrievsky et al. ( 1201 II) . To summarize, the temperature of the 
giants were inferred from their colours by adopting the calibra- 
tion of Alonso et al. (1 19961 1. and the temperature of the turnofF 
stars from the wings of the H„ line. The random error in the 
temperature is about 80K (Icr). For a given stellar temperature, 
the ionization equilibrium provides an estimate of the stellar 
gravity with an internal accuracy of about 0. 1 dex in log g, and 
the microturbulence velocity Vt can be constrained to within 
0.2 kms The adopted parameters are given in Table[3] Table 
|4] lists the uncertainty in the sulfur abundance originating in 
the errors in Jeff > log g , and v,. Because gravity is determined 
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Table 3. Adopted atmospheric parameters and NLTE S abundances for our sample of EMP stars. The solar abundance of sulfur 
has been taken from Caffau et al.: loge(S)o = 7.16, the NLTE Mg abundance (column 6) is from Andrievsky et al. (120101 1. The 
error (column 8) represents the uncertainty in the spectral fitting. The last column gives the number of sulfur lines used in the 
determination 
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3D+ 


3D+ 














NLTE 


LTE 


err. 


NLTE 


3D 


NLTE 


NLTE 


NLTE 




Star name 


Tcff 




Vt 


[Fe/H] 


[Mg/H] 


loge(S) 


+ 


loge(S) 


corr. 


[S/H] 


[S/Fe] 


[S/Mg] 


N 


Turnoff stars 


BS 16968-061 


6040 


3.8 


1.5 


-3.05 


-2.46 


4.90 


0.20 


4.33 


0.14 


-2.69 


0.36 


-0.23 




BS 17570-063 


6240 


4.8 


0.5 


-2.92 


-2.52 


5.05 


0.20 


4.61 


0.12 


-2.43 


0.49 


0.09 




CS 22965-054 


6090 


3.8 


1.4 


-3.04 


-2.42 


4.74 


0.20 


4.21 


0.13 


-2.82 


0.22 


-0.40 




CS 29499-060 


6320 


4.0 


1.5 


-2.70 


-2.14 


4.90 


0.20 


4.50 


0.12 


-2.54 


0.16 


-0.40 




CS 29506-007 


6270 


4.0 


1.7 


-2.91 


-2.25 


5.03 


0.20 


4.62 


0.12 


-2.44 


0.47 


-0.19 




CS 29506-090 


6300 


4.3 


1.4 


-2.83 


-2.22 


4.98 


0.20 


4.56 


0.12 


-2.48 


0.35 


-0.26 




CS 30301-024 


6330 


4.0 


1.6 


-2.75 


-2.25 


4.98 


0.20 


4.60 


0.12 


-2.44 


0.31 


-0.19 




Giants 


HD 2796 


4950 


1.5 


2.1 


-2.47 


-1.84 


5.25 


0.10 


4.97 


0.07 


-2.12 


0.35 


-0.28 


3 


HD 122563 


4600 


1.1 


2.0 


-2.82 


-2.19 


4.87 


0.10 


4.70 


0.06 


-2.40 


0.42 


-0.21 


2 


HD 186478 


4700 


1.3 


2.0 


-2.59 


-1.86 


5.07 


0.08 


4.84 


0.06 


-2.26 


0.33 


-0.40 


2 


BD+17 3248 


5250 


1.4 


1.5 


-2.07 


-1.39 


5.57 


0.10 


5.11 


0.10 


-1.95 


0.12 


-0.56 


2 


BD-18 5550 


4750 


1.4 


1.8 


-3.06 


-2.44 


4.65 


0.12 


4.41 


0.07 


-2.68 


0.38 


-0.24 


1 


BS 17569-049 


4700 


1.2 


1.9 


-2.88 


-2.11 


4.70 


0.20 


4.40 


0.06 


-2.70 


0.18 


-0.59 


2 


CS 22169-035 


4700 


1.2 


2.2 


-3.04 


-2.66 


4.50 


0.11 


4.20 


0.05 


-2.91 


0.13 


-0.25 


1 


CS 22186-025 


4900 


1.5 


2.0 


-3.00 


-2.39 


4.70 


0.10 


4.43 


0.06 


-2.67 


0.33 


-0.28 


2 


CS 22189-009 


4900 


1.7 


1.9 


-3.49 


-3.11 


4.25 


0.20 


3.97 


0.07 


-3.12 


0.37 


-0.01 


2 


CS 22873-055 


4550 


0.7 


2.2 


-2.99 


-2.34 


4.72 


0.10 


4.52 


0.10 


-2.54 


0.45 


-0.20 


1 


CS 22873-166 


4550 


0.9 


2.1 


-2.97 


-2.14 


4.70 


0.10 


4.55 


0.07 


-2.54 


0.43 


-0.40 


3 


CS 22878-101 


4800 


1.3 


2.0 


-3.25 


-2.56 


4.55 


0.16 


4.36 


0.05 


-2.75 


0.50 


-0.19 


1 


CS 22891-209 


4700 


1.0 


2.1 


-3.29 


-2.69 


4.40 


0.15 


4.21 


0.06 


-2.89 


0.40 


-0.20 


1 


CS 22892-052* 


4850 


1.6 


1.9 


-3.03 


-2.52 


4.48 


0.10 


4.29 


0.05 


-2.82 


0.21 


-0.30 


1 


CS 22896-154 


5250 


2.7 


1.2 


-2.69 


-2.06 


5.00 


0.08 


4.56 


0.07 


-2.53 


0.16 


-0.47 


2 


CS 22897-008 


4900 


1.7 


2.0 


-3.41 


-2.83 


4.35 


0.18 


4.10 


0.06 


-3.00 


0.41 


-0.17 


1 


CS 22948-066 


5100 


1.8 


2.0 


-3.14 


-2.59 


4.65 


0.14 


4.26 


0.06 


-2.84 


0.30 


-0.25 


1 


CS 22949-037* 


4900 


1.5 


1.8 


-3.97 


-2.42 


< 3.9 




< 3.6 


0.07 


< 3.7 


<0.5 


< -1.1 


2 


CS 22953-003 


5100 


2.3 


1.7 


-2.84 


-2.34 


4.72 


0.10 


4.41 


0.06 


-2.69 


0.15 


-0.35 


2 


CS 22956-050 


4900 


1.7 


1.8 


-3.33 


-2.59 


4.45 


0.20 


4.18 


0.06 


-2.92 


0.41 


-0.33 


1 


CS 22966-057 


5300 


2.2 


1.4 


-2.62 


-1.89 


5.12 


0.10 


4.62 


0.08 


-2.46 


0.16 


-0.57 


3 


CS 29491-053 


4700 


1.3 


2.0 


-3.04 


-2.34 


4.60 


0.10 


4.45 


0.05 


-2.66 


0.38 


-0.32 


3 


CS 29502-042 


5100 


2.5 


1.5 


-3.19 


-2.52 


4.50 


0.14 


4.13 


0.06 


-2.97 


0.22 


-0.45 


1 


CS 29518-051 


5200 


2.6 


1.4 


-2.69 


-2.02 


4.95 


0.10 


4.50 


0.07 


-2.59 


0.10 


-0.57 


2 


CS 30325-094 


4950 


2.0 


1.5 


-3.30 


-2.54 


4.64 


0.18 


4.19 


0.06 


-2.91 


0.39 


-0.37 


2 


CS 31082-001 


4825 


1.5 


1.8 


-2.91 


-2.09 


4.84 


0.10 


4.54 


0.07 


-2.55 


0.36 


-0.46 


3 


CS 22949-048 


4800 


1.5 


2.0 


-3.25 


-2.68 


4.42 


0.15 


4.07 


0.05 


-3.04 


0.21 


-0.39 


2 



* An asterisk after the name of the star means that the star is carbon-rich. 



from the ionization equilibrium, a variation in T^g will change 
log^ and also sometimes slightly influences v,. Hence, the to- 
tal error budget is not the quadratic sum of the various sources 
of uncertainties, but contains significant covariance terms (see 
Cayrel et al. 120041 for more details). 

Our \ogg value may be affected by NLTE effects (overi- 
onization) and the uncertainties in the oscillator strengths of 
the Fe and Ti lines. This was discussed in Cayrel et al. (120041 1. 
their section 3.1. : even for giants the spectroscopic \ogg and 
the logg value deduced from isochrones should not differ by 
more than about 0.3 dex. 



For the non-LTE computations of the sulfur abundance, 
Kurucz models without overshooting were used (Castelli et al., 
I1997I I. These models have been shown to provide LTE abun- 
dances very similar (within 0.05 dex) to those of the MARCS 
models used by Cayrel et al. (I20041 l and Bonifacio et al. (I2009I I. 

It has been shown (see Takeda et al. 120051 Korotin I2008I I 
that Mult. 1 of S I is affected by a significant negative non- 
LTE correction. To compute NLTE profiles of the S I lines, we 
used the model atom of Korotin (I2008I I. which contains 64 lev- 
els of S I and the ground level of S II. The radiative photoion- 
ization rates of all the levels taken into account are based on 
the new detailed ionization cross-sections listed in the Opacity 
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Table 4. Abundance uncertainties linked to stellar parameters. 
Case of a typical dwarf CS 29506-090, and a typical giant 
HD 122563. 



CS 29506-090 

A: refF=6300K, 
B: reff=6300K, 
C: reff=6300K, 
D: reff=6200K, 
E: reii=6200K, 



logg=4.3 dex, vt=1.4km s"' 
logg=4.2 dex, vt=1.4 km s"' 
logg=4.3 dex, vt=1.2km s"' 
logg=4.3 dex, vt=1.4km s"' 
logg=4.1 dex, vt=1.3 km s"' 



El. 




Ac-A 


Ao-A 


Af-A 


[Fe/H] 


-0.03 


0.03 


-0.06 


-0.07 


[S I/Fe] 


0.00 


-0.03 


0.09 


0.04 


[Fe I/Fe] 


0.02 


0.01 


-0.04 


0.01 


[Fe II/Fe] 


-0.01 


-0.01 


0.04 


-0.01 



HD 122563 

A: reff=4600K, 
B: r,ff=4600K, 
C: refr=4600K, 
D: reff=4500K, 
E: reff=4500K, 



log g= 1.0 dex, vt=2.0 km s"' 
logg=0.9 dex, vt=2.0 km s"' 
log g= 1.0 dex, vt=1.8 km s"' 
log g= 1.0 dex, vt=2.0 km s"' 
logg=0.6 dex, vt=1.8 km s"' 



El. 


Afi-A 


Ac-A 


Ad-a 


Af-A 


[Fe/H] 


-0.00 


0.06 


-0.09 


-0.06 


[S I/Fe] 


-0.02 


-0.02 


0.12 


0.04 


[Fe I/Fe] 


0.03 


0.03 


-0.11 


0.03 


[Fe II/Fe] 


-0.03 


-0.04 


0.11 


-0.03 



Project TopBase (hereafter TopBasefl The oscillator strengths 
of the lines are also taken from TopBase. The line profiles 
were computed from a modified version of the program MULTI 
(Carlsson 119861 Korotin et al. 1999). This version includes, 
in particular, opacity sources from ATLAS9 (Kurucz 1 9921) . 
which modify the continuum opacity and the intensity distribu- 
tion in the UV region. These sources are very important for de- 
termining of the radiative rates of the bound-bound transitions 
in the sulfur atom. In Fig. [3] we present the new non-LTE cor- 
rection computed for metal-poor stars for the 921.2 nm sulfur 
line of the Mult. 1 . This correction is sometimes rather differ- 
ent from the correction computed by Takeda et al. (2005). For 
example, for turnofF stars with a metallicity of about -3.0 and 
a temperature of 6000 K, the correction is found to be close to 
-0.6 dex, whereas it was only about -0.2 dex following Takeda 
et al. (I2005l l. 

We also computed the 3D correction (as explained in Caffau 
et al. 1201 Obi or Caffau & Ludwig I2007I I. This correction is 
small, but positive (Table O, and has been computed and ap- 
plied separately. 



4.1. Test of coherence 

With the model atom adopted here, there is good agreement 
(to within 0.1 dex) in the Sun and in Procyon between the 
S abundances derived from the different multiplets when 3D 
models are used. For Procyon, we found e(S) - 7.23 + 0.03 
from the 3D, LTE computations of Caffau et al. (I2007I I with 



a Korotin's NLTE correction (private communication). For the 
Sun, the mean sulfur abundance derived from these computa- 
tions is log e(S) = 7.16 + 0.05, for log e(H) = 12.0 (Caffau 
et al. (2010b). This value is very close to the meteoritic value 
(7.17 + 0.02, following Lodders et al.. |2009T l. The solar abun- 
dance loge(S) - 7.16 has been adopted as a reference in our 
computations. 

Moreover, for the most metal-rich star of our sample, 
BDh-17 3248, the S abundance could be determined from the 
lines of the Mult. 1 and the main line of Mult. 6 (Fig.[T^). The 
agreement between the NLTE S abundances derived from mul- 
tiplets 6 and 1 is very good when the 3D correction is also in- 
cluded (-0.02 for Mult. 6, +0.10 for Mult. 1): e(S) = 5.28 + 0.2 
from Mult. 6 and e(S) = 5.22 + 0.1 from Mult. 1. 
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Fig. 3. NLTE correction for the 921.2nm sulfur line in metal- 
poor stars as a function of the temperature, the gravity, and 
the metallicity for [Fe/H] =-2, [Fe/H] =-3 (dashed lines) and 
[Fe/H]=^ (solid line). This correction has been computed for 
[S/H] = [Fe/H]. 



5. Results 

In Table [3] we present the result of the computations for our 
sample of stars. The S abundances were derived from fits of 
synthetic spectra (Fig.|2|. The sulfur abundances could not be 
determined for the most metal-poor stars ([Fe/H] < -3.6 for 
giants, [Fe/H] < -3.2 for dwarfs), because the lines were too 
weak to be measured reliably. However, an upper limit to the 
sulfur abundance was given for CS 22949-037 because this 
upper limit can help us to explain the chemical anomalies of 
this star (see section l53T l. The error in the sulfur abundance, 
given in Table [3] represents the uncertainty in the spectral fit- 
ting. This error strongly dominates the total error in both [S/Fe] 
and [S/Mg]. 

5.1. [S/Fe] vs. [Fe/H] and comparison to previous 
studies 

Fig. m shows [S/Fe] as a function of [Fe/H] for our sample of 
EMP stars. Below [Fe/H] =-2.9, the slope of the relation is very 
small : a - -0.012 + 0.006, thus [S/Fe] is practically constant 
with a mean of [S/Fe] = -1-0.35 + 0.10. The level of this plateau 
is in good agreement with the value found by Nissen et al. 
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Fig. 4. Final NLTE [S/Fe] ratios for our BMP stars. Solid sym- 
bols show the turnofF stars, open symbols the giants. 



(120041 12007] l and Caffau et al. (l2005T l. [S/Fe] 
for the most metal-poor stars of their sample). 
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Fig. 5. Spectrum of BD-18 5550 in the region of Mult. 1 of 
S I. Only the main line can be measured in our spectrum, as the 
others are severely blended with telluric lines. The thick line 
shows the theoretical profile computed with our value of the S 
abundance (Table 13). The profile computed with the abundance 
adopted by Takeda & Takada-Hidai (2010) (thin line), is not 
compatible with the observed spectrum. 



• Behavior of stars with -2.5 < [Fe/H] 
It is interesting that none of the stars in Fig.|4]is found to be "S- 
rich" ([S/Fe] > +0.7), in contrast to several stars in the interval 
-2 < [Fe/H] < -1 (IsraeUan & Rebolo[200ll Takada-Hidai et 
al. l2052l Takeda et al. 12005. and Caffau et al. l2005T l. 

The most S-rich star in Israelian & Rebolo ( 1200 11 1 is 
HD 2796, for which [S/Fe]=-i-0.8. It is included in our sample, 
but its very high S abundance is not confirmed by our measure- 
ment of the lines of Mult. 1 from which we find [S/Fe]=H-0.3. 
Israelian & Rebolo found that EW=0.26pm for the main line 
of Mult. 6 in this star, but although our spectrum has a very high 
S/N in this region (over 500/pix at 869 nm), we do not detect 
the line (Fig. [TJ)). Its equivalent width is clearly smaller than 
0.26 pm, the error in the equivalent width was probably un- 
derestimated by Israelian & Rebolo. Ryde & Lambert (I2004I I. 
reanalysed three other S-rich stars from the list of Israelian & 
Rebolo with 0.55 < [S /Fe] < 0.7 and for all of them they found 
also found a normal sulfur abundance. 



In Caffau et al. (120051) . three other stars were found to be S- 
rich (with [S/Fe] > -hO.7), BDh-02263, BD-h05 3640 and CD- 
45 13178 (HD 181743), their metallicities being in the interval 
-2.2 < [Fe/H] < -1.1, thus outside the range of metallicity 
explored in our paper Moreover, in HD 181743, Nissen et al. 
(120071) found from the Mult. 1 a normal abundance of sulfur 
([S/Fe]=+0.34), but on the other hand, in BD+05 3640 the high 
sulfur abundance has been confirmed (but at a lower level) by 
Caffau et al. (I2010al l by the analysis of the Mult. 3. 

It remains unclear wether S-rich stars exist with metallici- 
ties in the interval -2.5 < [Fe/H] < -1.1. 

• Behavior of stars with [Fe/H] < -2.5 
Recently, Takeda & Takada-Hidai (120101 1 measured the Mult. 3 
of S I in metal-poor stars and found [S/Fe] to be constant in the 
interval -2.5 < [Fe/H] < -1.1, but to increase suddenly for 
[Fe/H] < -2.6. They have only five stars with [Fe/H] < -2.6, 
and in one of them the S line cannot be detected. For G64-37, 
the uncertainty is so high that the detection is in doubt (see 
their Fig. 13). Of the three remaining stars, two of them are in 
common with our sample: HD 122563 and BD-18 5550. For 
these stars, we find "normal" S abundances of [S/Fe] =-1-0.42 
and H-0.38 from Mult. 1 (TableO. 

To understand the cause of this discrepancy, we remeasured 
the equivalent widths of the 1045.5449 nm line in the spectra of 
Takeda & Takada-Hidai (2010), neglecting the other two lines, 
which are even fainter. The lines of Mult. 3 used by Takeda & 
Takada-Hidai are weaker than those of Mult. 1 and their spec- 
tra have only half the resolution of ours. For HD 122563, we 
find that EW=1.5 pm, corresponding to [S/Fe]= H-0.45, in ex- 
cellent agreement with our determination. Since the sulfur lines 
are weak, the value of [S/Fe] found by Takeda & Takada-Hidai 
was simply corrected by A[S/Fe] - log(EW2/EWi), where 
EWl and EW2 represent the old and new equivalent widths. 
BD-18 5550 is the most S-rich star in the sample of Takeda 
& Takada-Hidai with [S/Fe] =-1-0.76, against our value of +0.38 
from Mult. 1 (see Fig. |5]l. For this star, the S/N of the spec- 
trum of Takeda & Takada-Hidai seems to be rather high, but 
we note that the full width at half maximum (FWHM) of the 
lines is much larger in this spectrum than in the spectrum of 
HD 122563, while in our UVES spectra the FWHM of the lines 
in both stars is about the same. If we fix the FWHM of the lines 
in BD-18 5550 to the same value as in HD 122563, we find the 
equivalent width of the 1045.5449 nm line to be about 1.2 pm 
and [S/Fe] =+0.62. Adopting the Cayrel formula (Cavrel [T9g8] l 
to estimate the precision of this measurement (with S/N =150), 
we find the error of this equivalent width to be close to 0.5 pm. 
Hence, the error in the abundance ratio should be about 0.2 dex, 
which makes their result compatible with ours. 



5.2. [S/Mg] vs. predicted supernova yields 

It is often noted that magnesium should be a more reliable ref- 
erence element than iron, because Mg is only formed in mas- 
sive SN II and its production is less affected by explosive nu- 
cleosynthesis, mixing, and "fallback". As a consequence, the 
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Fig. 6. [S/Mg] vs. [Fe/H] (top) and [Mg/H] (bottom). Symbols 
as in Fig.m 



detail when we have completed a new precise analysis of disk 
stars, including NLTE and 3D corrections. 
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Fig. 7. [X/H] abundances in CS 22949-037 against the solar 
condensation temperatures (Lodders, I2003I I. There is no cor- 
relation between Tcond and the abundance of the elements, in 
particular S and Zn have a condensation temperature lower than 
Na and they are more deficient, at variance with the abundances 
in A Boo stars. (The filled circles represent the abundances cor- 
rected or not sensitive to NLTE effects, the open symbols are 
abundances that could not be corrected for NLTE effects) 



predictions for [S/Mg] should be more robust than for [S/Fe] 
(Woosley & Weaver. [T9951 Shigeyama & Tsujimoto 1998J . 

The NLTE abundance of Mg for our sample of stars was 
computed by Andrievsky et al. (2010). Fig. |6] shows the ob- 
served value of [S/Mg] as a function of [Fe/H] and [Mg/H]. The 
scatter is small, but two stars seem to have a high [S/Mg] ratio: 
BS 17570-063 ([Fe/H]=-2.92) and CS 22189-009 ([Fe/H]=- 
3.49). In both cases, the high [S/Mg] ratio seems to be real, 
and not the consequence of an uncertainty in the magnesium 
abundance (see Fig. |6]). The full spread of NLTE [Mg/Fe] for 
the metal-poor stars below [Fe/H] < -1 is definitely larger 
than expected from the error estimate of the analysis itself (see 
Gehren et al.. 120061 Andrievsky et al., 2010, and also Cayrel et 
al.. l2004l l and the anomalies of the magnesium abundances are 
not clearly correlated with the other elements. 

In the interval -3.4 < [Fe/H] < -2.5, we can define a mean 
value of [S/Mg] as done for O, Na, Al and K by Andrievsky et 
al. ( [20T0I 1: [S/Mg] = -0.32 + 0.14. In the same metallicity 
interval, the scatter in [S/Fe] around the mean is only 0.12 dex 
(see FigHJl. It appears that, as already observed in LTE (Cayrel 
et al. l2004l l. the scatter of the abundance ratios is slightly larger 
when Mg rather than Fe is used as the reference element. This 
is surprising because Fe forms in processes quite different from 
those forming Mg and S (explosive burning for Fe, and mainly 
hydrostatic C, Ne, and O burning for Mg and S). 

For [Mg/H] < -2.7, a slight increase in [Na/Mg], [Al/Mg], 
and [K/Mg] as [Mg/H] decreases was seen by Andrievsky et al. 
(1201 Oi l. A slight increase in [S/Mg] is also possible (Fig.|6]l, but 
there are very few stars with [Mg/H] < -2.7 in our sample. We 
plan to discuss the evolution of [S/Mg] in the Galaxy in more 



5.3. The case of the peculiar C-rich star 
CS 22949-037 

CS 22949-037 is a very peculiar giant: it is extremely metal- 
poor ([Fe/H] ^ -4.0), and its unusual chemical composition 
was first pointed out by Mc William et al. ( i 1995 J and later stud- 
ied in more detail by Norris et al. ( 1200 11 1 and Depagne et al. 
(2002). This star exhibits light element (C-N-O, Na, Mg) en- 
hancements that are quite similar to those associated with the 
hyper iron-poor stars HE 0107-5240 and HE 1327-2326 with 
[Fe/H] < -5 (see Frebel et al., l2005l l2008b . It is generally 
admitted that these stars are the product of specific zero-metal 
supernovae. However, Venn & Lambert (120081 ) argued that this 
abundance pattern might be associated with the same peculiar- 
ity that is used to account for ABoo stars (presumably related 
to dust depletion). If this were so, the abundance of the el- 
ements should be smoothly correlated with the condensation 
temperature Tcond- The abundances of S and Zn could help us 
to choose between the two hypotheses, because their temper- 
atures are intermediate between those of O and Na, both of 
which are strongly enhanced in CS 22949-037. If dust deple- 
tion is responsible for the pattern of CS 22949-037 we would 
expect that sulfur and zinc also be strongly enhanced in this 
star. The abundance of zinc was measured in CS 22949-037 
by Depagne et al. (120021 ). The sulfur lines are not visible in 
this star but (as can be seen in Table [3]) an upper limit can be 
estimated as [S/H] < -3.5 or [S/Fe] < H-0.5. In Fig. |7] we 
have plotted [X/H] vs. Tcond- The abundances [X/H] of C, N, 
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O, Ca, Fe, and Zn are taken from Depagne et al. (I2002I I: for 
the other elements Na, Mg, Al, S, Sr, and Ba, the NLTE abun- 
dances (Andrievsky et al.. . 20071 l2008ll2009ll20T0ll20TTi) have 
been preferred. The oxygen abundance, determined from the 
forbidden line, is insensitive to the NLTE effects. It can be seen 
that CS 22949-037 does not show any clear signature of dust 
depletion. 

5.4. Comparison with predicted SN yields 



Table 5. [Zn/H] and [S/Zn] in the most metal-poor DLAs. 
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Fig. 8. Comparison of the most recent O, Na Mg, Al, S, and K 
abundances in BMP stars to the predicted yields of metal-poor 
SNe and hypernovae. The solid lines represent the predictions 
of different models, and the dots the mean values of the ob- 
served abundances in the interval -3.0 < [Mg/H] < -2.0 



In Fig. [8] we compare the abundance ratios of the light 
metals to some predicted abundances in the ejecta of metal- 
poor SNe or hypernovae following Woosley & Weaver (|19951 l. 
Kobayashi et al. ('2006^, and Heger & Woosley (|20T0] |. The 
mean values of the observed abundances of the different ele- 
ments are computed in the interval -3.0 < [Mg/H] < -2.0, 
their uncertainty being of the order of 0. 1 dex. The amount of 
ejected S seems to be generally overestimated by the models. 
If we consider all the elements from O to K, the best agree- 
ment is obtained with the predictions by the models of Heger 
& Wooslev l2010l (hvpotheses B or D, with a low mass cutoff at 
IOMq), mainly because their predictions show a better agree- 
ment with the observed K abundance. 

5.5. [S/Zn] in metal-poor stars and damped Lya 
systems 



o o o ^ 
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Fig. 9. [S/Zn] vs. [Zn/H] in BMP stars (circles) and in metal- 
poor DLAs (triangles). 
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Nissen et al. ( I2004I I suggested that the [S/Zn] ratio is higher 
in stars than in damped Lyo- systems. Neither S nor Zn shows 
much affinity for dust, although MgS has been recently pro- 
posed as a possible grain constituent by Zhang et al. ( I2009I I. 
Thus, depletion onto dust in DLAs probably cannot explain this 
difference. Nissen et al. (2007) have shown that in the interval 
-2.5 < [Zn/H] < this effect is considerably reduced if one 
adopts the non-LTB corrections to the abundances of S and Zn. 
Fig. |9] shows [S/Zn] vs. [Zn/H] for our sample of metal-poor 
stars and for the most metal-poor DLAs with both Zn and S 
abundance determinations (Table |5]l. Since the Zn line is very 
weak in turnoff stars, we could only measure Zn abundances 
in our giant stars (Cayrel et al.. l2004l i. Following Takeda et al. 
(12005 non-LTB effects for Zn are not very important in these 
stars (< +0.1dex), but we have applied the small correction ob- 
tained by interpolation in their tables. 

With stellar NLTB effects taken into account for both S and 
Zn, Fig. |9] shows that the sulfur abundance in metal-poor stars 
is rather similar to that in DLAs - if anything slightly lower. 
However, the chemical evolution of the Galactic halo might still 
be different from the average of DLAs. 

In BMP stars as in DLAs, the [S/Zn] ratio is approximately 
solar. This is due to positive values of both [S/Fe] (Fig.|4]i and 
[Zn/Fe] at very low metallicity (Cayrel et al. .2004. see also 
Primas et al. 120001 and Nissen et al. l2004l l. We will discuss the 
evolution of the [S/Zn] ratio in the Galaxy in more detail when 
we have completed a new precise analysis of S abundances in 
less metal-poor stars. 



6. Conclusions 

We have determined S abundances in 33 metal-poor stars, in- 
cluding 21 with [Fe/H] < -2.9, and find [S/Fe] to be fairly 
constant at [S/Fe] - 0.35 + 0.10. However, a slight increase 
in [S/Mg] when [Mg/H] decreases cannot be excluded (Fig.|6]l. 
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The same effect is also observed for Na, Al, and K (Spite et al. 
120091 Andrievsky et al. l20T0] i. 

We do not confirm earlier claims of a significant rise in 
[S/Fe] with decreasing [Fe/H], but corroborate the findings of 
Nissen et al. (l2004l |2057] | and Ryde & Lambert (l2004b that S 
behaves like the other a-elements, also at very low metallicity. 

The low sulfur abundance (along with a low abundance 
of zinc) in CS 22949-037 suggests that the abundance pat- 
tern in this star (and thus by analogy in the hyper iron-poor 
stars HE 0107-5240 and HE 1327-2326) is the result of spe- 
cific zero-metal supernovae and is unrelated to dust depletion 
as suggested by Venn & Lambert (I2008I I. 

When the NLTE effects are taken into account, the pat- 
tern of the light elements from O to K relative to magnesium, 
is rather well represented by ejecta of zero-metal supernovae 
(Heger & Woosley|20l0l). 

The ratio [S/Zn] is found to be approximately solar as it is 
in DLAs since [S/Fe] and [Zn/Fe] are both equally positive at 
low metallicity. 
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